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. Abstract. The Color Glass Condensate formalism and its application to high energy 

I heavy ion collisions at RHIC are discussed. We argue that the RHIC data supports 

the view that the Color Glass Condensate provides the initial conditions for gold- 
gold collisions at RHIC while final state (Quark Gluon Plasma) effects are responsible 
' for the high pt suppression in mid rapidity. At forward rapidities in deuteron-gold 

collisions, however. Color Glass Condensate is the underlying physics of the observed 
suppression of the particle spectra and their centrality dependence. 
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^ ' 1. Introduction 



The recent data from forward rapidity deuteron-gold collisions at the Relativistic 
Heavy Ion Collider (RHIC) in Brookhaven National Laboratory (BNL) has generated 
^ ! much interest and has led to intense debate on the underlying physics of the forward 

^ ' rapidity processes in deuteron-gold collisions. While the conventional models of multiple 



X 



scattering based on Glauber theory predicted an enhancement of the Cronin effect in 
the forward rapidity region, the data shows a clear suppression of the negatively charged 
hadron spectra, as compared to the normalized proton-proton collisions, in agreement 
with the predictions based on the Color Glass Condensate formalism. The most recent 
data on centrality dependence of the suppression in the forward rapidity region also 
agrees with the predictions of the Color Glass Condensate formalism and is in clear 
contradiction with expectations based on the conventional models. 

The Color Glass Condensate formalism is an effective theory of QCD at high energy 
(small x). At high energy or equivalently (for fixed Q'^) at small x, gluons are the most 
abundant partons. While one can not calculate parton distribution functions in QCD 
using perturbation theory, one can calculate their evolution with x or in pQCD. 
The evolution of the gluon distribution function is described by the DGLAP evolution 
equation which sums log type corrections to the gluon distribution function. In 
DGLAP formalism, the gluon distribution function xG{x, Q^) grows fast with both x 
and as shown in Eq. for fixed a coupling constant a^. It is worth noting that 
DGLAP evolution equation is just the renormalization group equation for the parton 
number operator. 
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xG{x, Q'^) ~ exp 



as logl/xlogQ^ 



(1) 



If the phase space in is hmited and there is a large phase space available in 
then BFKL is the applicable formalism. In BFKL approach, one sums terms of the form 
Us log 1/x which leads to an even faster grows of the gluon distribution function with x 
(energy) compared to the DGLAP gluons 



xG'(x,g2) ^ e^'°e^/^ (2) 

with A = ^^^A log 2 for a fixed coupling constant. In other words, pQCD radiation leads 
to a fast growth of the gluon distribution function with energy or x which makes a hadron 
a dense system of gluons at high energy. The gluon density in a high energy nucleus is 
even larger by a factor of ~ A^^^ due to the Lorenz contraction of the longitudinal size 
of the nucleus at high energy. However, this fast growth can not go on forever since it 
would lead to violation of unitarity for physical cross sections. 

Both DGLAP and BFKL approaches arc linear in the gluon density. In other 
words, they only include radiation of gluons. As the phase space density of gluons 
increases, it becomes as probable for a gluon at small x to recombine with another 
gluon at small x and become a gluon at a higher x as it is for a gluon at larger x to 
radiate another gluon at a smaller x. This recombination of gluons slows down the 
growth of the gluon distribution function and leads to its eventual unitarization. This 
was first investigated by Gribov, Levin and Ryskin (GLR) and later on by Mueller and 
Qiu (MQ). The evolution equation for the gluon distribution function is written as (in 
the Double Logarithm Approximation) 

xG{x,Q )-7775 — --( ) -^x G {x,Q )(3) 



9 log 1/x (9 log g2 ' ' TT 'Q^ 

Here G'^{x, Q^) is the four point function of gluons while the standard gluon distribution 
function xG{x, Q'^) is the two point function of gluons. It is quite customary to 
factorize the four point function G^ in terms of the gluon distribution function G as 
G^(a;, Q'^) ~ ;^[G(a;, Q^)]^ where i? is a scale which can not be determined from pQCD. 
This factorization is an approximation where one ignore gluon correlations, valid at large 
Nc- The fact that this is a higher twist correction to the gluon distribution function is 
manifested in the presence of the factor in the non linear correction. 

As the phase space density of gluons becomes large, even the GLR-MQ approach 
breaks down since all higher point functions become as large as the two point function 
and one must keep all higher order terms in the evolution equation. The Color Glass 
Condensate is an effective field theory approach to QCD at high energy (small x) which 
extends the applicability of pQCD to a dense system of gluons. 
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2. The Color Glass Condensate 

The effective action for QCD at small x is given by (for a recent review, see PP) 
S = tJ d\tF[p%Xt)] 

- J d^x^G^ + d^Xtdx~5ix-)p''ixt)trTaW.^,^[A-]{x-,Xt) (4) 
where G^"^ is the gluon field strength tensor 

Gr = o'^a: - d'^A': + gfa,cA^,A: (s) 

Ta are the SU{N) color matrics in the adjoint representation, and W is the path ordered 
exponential along the x^ direction in the adjoint representation of the SU{Nc) group 



W-oo,oo[A ]{x ,Xt)=Pexp 



ig / dx'^A^ {x , Xt)Ta 



(6) 



In order to calculate a physical quantity, one averages over the gluonic field configuration 

I[Dp''][DAii]0{A) exp{zS[p,A]} 

J[Dp-][DA^a]exp{zS[p,A]} ^' ^ 

The weight function 

Z = exp{- J d^XtFlp"]} (8) 

appearing in ((Tj) is the statistical weight of a particular configuration of the two 
dimensional color charge density p°'{xt) inside the hadron. Initially, this statistical 
weight is taken to be a Gaussian 

FWt)] = ^Pli^t) (9) 

The Gaussian approximation is valid as long as the color charge density is large and 
random. Also, if one is interested in the gluon correlation at the same rapidity, one 
needs to include the longitudinal structure of the sources into account. This is done in 
[2] where one solves the classical equations of motion in the presence of sources which 
have a longitudinal size. The result for the two point function is 

Gy{y, xr, y', xl) = {A{y, x^)A,{y', x'^)) (10) 

Performing the color averaging with the extended sources, we get 



^ 4(iV2 - 1) 



N,x 



2 



(^I^qcd) 



(11) 



The momentum space distribution of the gluon distribution function obtained from the 
Fourier transform of Eq. (jllj) has a much milder divergence at low kt than the pQCD 
distributions due to the non linearities of the gluons included in a classical description. 
One can also show that gluons typically have a transverse momentum Qs{x, bt,A), called 
the saturation scale, which can be much larger than Aqcd- This is shown in Fig. 
which illustrates the fact that most gluons reside in a state with a momentum ~ Qs 
described by a strong classical field. As one goes to smaller x, one needs to include 
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Figure 1. Phase space distribution of gluons 



quantum evolution effects, namely, sum the a^logl/a; terms which become large. This 
has been done and leads to a non linear functional equation for the weight functional 
F{p) which can then be used to write an evolution equation for any desired physical 
observable. The JIMWLK equation for the weight functional [3] can be written in a 
compact form as (first written in this form by H. Weigert) 

dyZy[U] = -\ I VI xtn ' KMU] (12) 

with y = logl/x. These equations are highly non linear and coupled so that one has 
to develop approximate techniques to solve them. They have also been investigated 
numerically. One particular limit which simplifies the non linear JIMWLK equations, is 
the large limit which was used in to derive an equation for the dipole scattering 
cross section. This equation is known as the BK equation and is widely regarded as the 
simplest possible non linear equation which respects perturbative unitarity. The BK 
equation is simplest in the momentum space where it can be written as 

= 2axi~jJ-^)Nix, h) - aN'ix, fc„ h) (13) 
oy o log k 

where N is related to the dipole cross section via 

iV(x, ku bt) = f p\e-'''-'-^aix, n, h) (14) 

and the dipole cross section a{x,rt,bt) is defined in terms of the expectation value of 
the correlator of two Wilson lines 

a{x,n,bt) = j^Tr < 1 - U{xt)U\yt) > (15) 

with Tt = Xt — yt and bt = {xt + yt)/'^- The BK equation has been investigated both 
numerically and analytically and approximate analytic solutions as well as complete 
numerical solutions have become available. 
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3. Applications to RHIC 

The Color Glass Condensate was first applied at RHIC to the measured hadron 
multiplicities in nucleus-nucleus collisions at -y/i = 130 GeV The starting point 
is the kt factorized form of gluon production cross section in heavy ion collisions 

da AttNc 1 /■ ,,2 / 7 2\ / / 7 \2\ 

d^ = j\j2 J "^t )(I^{X2, [Pt - kt) ) (16) 

where Xi,X2 = -^exp(±?7). To get particle multiplicities, one integrates (fTB|) over 
the transverse momentum and divides by the inelastic cross section of nucleus-nucleus 
collisions. The result is shown in Fig. (j2I) where a good agreement with the data at 
^/s = 130 GeV is obtained. The centrality and energy dependence are also in good 
agreement with the RHIC data. It should be noted that multiplicities are dominated 
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Figure 2. Hadron multiplicity 

by low transverse momenta. For example, more than 95% of produced particles are 
produced with transverse momentum less than 1 GeV which indicates that multiplicities 
probe xi,X2 of the nucleus wave functions which are smaller than 0.01 at mid rapidity. 

Eq. ()16p was also applied to hadron transverse momentum spectrum in mid rapidity 
nucleus-nucleus collisions at ^/s = 130 GeV and a reasonable agreement was obtained. 
This would indicate that the observed suppression of the hadron spectra in mid rapidity 
nucleus-nucleus collisions is due to initial state (and not Quark Gluon Plasma) effects. 
Deuteron-gold collisions were performed in order to distinguish initial state from final 
state effects. The observed back to back hadron correlations and a lack of suppression 
in hadron spectra in mid rapidity deuteron-gold collisions basically rule out high gluon 
density effects at high pt in mid rapidity gold-gold collisions and are a strong indication of 
formation of a dense, strongly interacting medium which is most likely partonic in origin. 
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As such, the Color Glass Condensate provides the initial conditions for the formation 
of the Quark Gluon Plasma P . However, the physics of mid rapidity nucleus-nucleus 
collisions is mostly final state interactions between produced partons. Therefore, it is 
not the best place to look for manifestation of the Color Glass Condensate because, at 
high pt, the final state effects in mid rapidity nucleus- nucleus collisions are clearly more 
important than initial state effects. 

On the other hand, the forward rapidity region in deuteron-gold collisions is mostly 
free of these final state interactions since there is no Quark Gluon Plasma expected to be 
formed there. Also, in the forward rapidity region in deuteron-gold collisions, one probes 
the smallest x in the nucleus wave function kinematically possible. This is where the 
Color Glass Condensate effects are expected to be the strongest. Another reason why the 
forward rapidity region is the best place to look for the Color Glass Condensate effects in 
a collider environment is that in forward rapidity, the experimental pt coverage is limited 
by the kinematics. Therefore, evolution in pt is not very important while the effects of 
X evolution are maximized. This is crucial since the Color Glass Condensate formalism 
does not include a^logQ^ effects (summed by DGLAP) which become important at 
high Pt (the double log region where one sums as log log 1/x terms is indeed included 



In the forward rapidity region of deuteron-gold collisions one probes the large x 
part of the deuteron wave function. To simplify the discussion, we will focus on proton- 
nucleus collisions. In this kinematic region, proton is a dilute system of partons with 
their distributions well described by pQCD. On the other hand, the nucleus is probed at 
very small x so therefore, we treat it as a Color Glass Condensate. The starting point 
is to consider scattering of a valence quark of the proton on the nucleus |2j. One could 
also consider radiation of a gluon from the quark scattering on the nucleus [H]. The 
properties of these cross sections have been studied in detail in CD] where it shown 
that Color Glass Condensate formalism, at the classical level, leads to the Cronin effect 
while inclusion of x^j evolution leads to a suppression of the nuclear enhancement factor. 
Since the case when the valence quark scatters on the nucleus is simpler (and possibly 
the most important in the forward rapidity and low pt), we consider scattering of a quark 
on the nucleus described as a Color Glass Condensate. The differential cross section is 
given by 



where the dipole cross section a{x,rt,bt) is given by (fTK|l . It is interesting to note that 
this cross section is finite as unlike the pQCD cross sections which are divergent. 
To obtain cross sections for hadron production, one needs to convolute ()17|1 with the 
quark distribution function in a proton and a quark-hadron fragmentation function 



In Fig. we show the result of a recent calculation JT] based on ()18|) for minimum 
biased deuteron-gold collisions in forward rapidity {y = 3.2) and compare it with the 



in CGC). 




(17) 




(18) 
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nuclear modification factor RaA for negatively charged hadrons measured by BRAHMS 
[T^ (our calculation is for the sum of positively and negatively charged hadrons, but 
in the kinematic region considered here, this does not make a big difference). The 
agreement at low pt is quite striking specially since there is absolutely no free parameter 
used here. 
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Figure 3. Minimum bias RdA for forward rapidity {y = 3.2) deuteron-gold collisions. 

The physics of particle production in forward rapidity and low pt can be understood 
in simple terms. The valence quark of the deuteron scatters from the fully developed 
(as allowed by the kinematics) wave function of the nucleus, characterized by Qf{y) = 
^e^yl'^ ~ 2.2 GeV (at y = 3.2) which is roughly how much transverse momentum the 
scattered quark accumulates. Assuming < z >~ 0.8 in the fragmentation process, this 
means quark-nucleus scattering can describe particle production up to ~ 1.8 GeV. As 
one goes to higher pt, one needs to include radiation of gluons in analogy with DGLAP 
radiation. In [Hj, gluon radiation from a quark scattering on the nucleus is calculated 
and must be included at higher pt. Indeed, the suppression of hadron multiplicities in 
deuteron-gold coUisions in the forward rapidity region had been predicted jTHI based 
on the behavior of the gluon production cross section in different pt regions as was the 
striking centrality dependence of the hadron spectra in different rapidity regions. The 
fact that the cross section for scattering of quarks on the nucleus has similar properties 
is not surprising since both quark scattering and gluon production cross section depend 
on the (fundamental vs. adjoint) dipole cross section. 

4. Summary 

The recent results from the Relativistic Heavy Ion Collider have provided much 
excitement and enthusiasm in the high energy heavy ion community. On the one hand, 
there is overwhelming evidence for the formation of a deconfined state of matter, most 
likely of partonic nature in thermal equilibrium, in mid rapidity heavy ion collisions. On 
the other hand, we have, very likely, seen the first unambiguous evidence for the Color 
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Glass Condensate in the forward rapidity deuteron-gold collisions. Even though the 
Color Glass Condensate correctly predicted rapidity, energy and centrality dependence 
of hadron multiplicities in gold-gold collisions, the more conventional models, based on 
Glauber multiple scattering, were also able to accommodate the data. 

The spectacular failure of the conventional models applied to the forward rapidity 
region in deuteron-gold collisions and the agreement of the Color Glass Condensate 
predictions with the observed suppression of the nuclear modification factor and its 
centrality dependence have provided strong evidence in favor of the Color Glass 
Condensate formalism. To put this on a more firm footing experimentally, one needs to 
measure more observables, such as photons, dileptons, etc. ^Sl^^- Another deuteron- 
gold run at RHIC will go a long way towards establishing the properties of the Color 
Glass Condensate and is urgently needed. 
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